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Abstract. The insulating state of magnetite (Fe 3 04) can be disrupted by a 
sufficiently large dc electric field. Pulsed measurements are used to examine the kinetics 
of this transition. Histograms of the switching voltage show a transition width that 
broadens as temperature is decreased, consistent with trends seen in other systems 
involving "unpinning" in the presence of disorder. The switching distributions are 
also modified by an external magnetic field on a scale comparable to that required to 
reorient the magnetization. 
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1. Introduction 

Magnetite is an archetypal strongly correlated transition metal oxide, with properties 
not well described by single-particle band structure. Below 858 K, magnetite, which 
may be written as Fe^ + (Fe 2+ Fe 3+ )B04, is ferrimagnetically ordered, with the A and 
B sublattices having oppositely directed magnetizations. The moments of the five 
unpaired 3d electrons of the tetrahedrally coordinated A-site Fe 3+ ions are compensated 
by those of the octahedrally coordinated B-site Fe 3+ ions. The net magnetization results 
from the octahedrally coordinated B-site Fe 2+ that have four unpaired 3d electrons [I]. 
Upon cooling, bulk magnetite undergoes a first-order phase transition from a moderately 
conducting high temperature state to a more insulating low temperature state at what is 
now called the Verwey[2] temperature, Ty ps 122 K. The change in electronic properties 
is coincident with a structural transition from a high temperature cubic inverse spinel 
to a low temperature monoclinic unit cell. The nature of the ordered insulating state 
remains an active topic of current research[3j HI El [6]. Experiments indicate the onset 
of multiferroicity[7] in magnetite below 40 K[8], further highlighting the rich physics in 
this correlated system. 

Recently, nanostructured electrodes have been used to apply strong electric fields 
in the plane of magnetite films[9] [10]. Below Ty, a sufficient applied voltage triggers a 
breakdown of the comparatively insulating low-temperature state and a sudden increase 
in conduction [9j [10] . This is an example of electric field-driven breakdown of a gapped 
state in strongly correlated oxides [TT1 [121 [T3] similar to Landau-Zener breakdown in 
classic semiconductors. The electric field-driven transition in magnetite is consistent 
with expectations [13] based on such a mechanism (via geometric scaling [9], [10], lack of 
intrinsic hysteresis [II] . changes of both contact and bulk resistance at the transition 
[T5| [T6]). These prior experiments examined films of various thicknesses, from 30 nm 
to 100 nm. No strong thickness dependence was observed in the switching properties, 
consistent with the applied lateral electric field at the sample surface acting as the driver 
of the breakdown (though thinner films showed a less pronounced Verwey transition in 
low-bias resistance vs. temperature measurements, consistent with expectations). 

Here we report studies of the statistical variations of this electric field-driven 
transition in Fe3C>4, as a function of temperature and magnetic field perpendicular 
to the film surface (out-of-plane) . We find that there is a statistical distribution of 
switching voltages, Vsw, that becomes more broad and shifts to higher voltages as T is 
reduced. We discuss these trends in the context of switching kinetics in other systems 
that exhibit similar trends. The application of a magnetic field perpendicular to the 
plane of the Fe3C>4 film alters the Vswj shifting the mean by several mV (several percent) 
and changing its shape, within a range of fields comparable to that required to reorient 
the magnetization out of plane. 
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2. Experimental Techniques 

The 50 nm Fes04 (100) thin films used in the present study were grown on (100) oriented 
MgO single crystal substrates as described elsewhere [171 [18]. Contact electrodes (2nm 
adhesion layer of Ti and 15 nm layer of Au) were patterned by e-beam lithography on the 
surface of the Fe3C>4 film. As before [HUH], Vsw scales linearly with the channel length, 
L (the electrode spacing), implying an electric field-driven transition. Long channels 
(L > 100 nm) required large switching voltages that would alter the electrode geometry 
over numerous switching cycles, distorting the shape of Vsw histograms. To minimize 
Vsw, electrodes separated by 10-30 nm were patterned using a self-aligned technique 
|10j . Electrical characterization of the devices was performed using a semiconductor 
parameter analyzer (HP 4155A). To minimize self- heating when in the conducting state, 
the voltage was applied as pulses 500 /is in duration with a 5 ms period [TU [15]. The 
samples were cooled below Ty with no magnetic field applied, and the distribution of 
Vsw was obtained by executing several thousand consecutive forward pulsed J-V sweeps 
in the vicinity of the transition point (typically a 0.2-0.3 V range) at a fixed (to within 
50 niK) temperature, and recording the number of switching events at each voltage. 

Each voltage value is essentially an independent test to see if switching takes place 
under the pulse conditions. Hence, some sweeps show one (figured] a) or several (figure [1] 
b) switching events. Even if the system is switched to the conducting state at Vsw(l), 
it may return to the Off state between pulses, and then switch to the On state at some 
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Figure 1. (color online) Details of Vsw distribution experiment, (a and b) The 
fragments of I-V curves in the vicinity of a transition demonstrating one (a) 
and three (b) switching events in a single pulsed I-V cycle, (c) Typical SEM 
image of Ti/Au electrodes, patterned on magnetite film surface, separated by 
nanogap < 100 nm. (d) An example of Vsw distribution histogram at 90 K. 
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higher voltage, Vsw (2), and so on. The Vsw distribution at a particular temperature 
is built by recording all switching events over several thousands (3000-6000) of I-V 
cycles and then counting the number of switchings at a certain Vsw, to produce a 
of counts" vs. Vsw histogram. A typical V"sw distribution at 90K is shown in figure [1] 
d. The distribution is a single peak, symmetrical around the most probable Vsw value. 

3. Results and Discussion 

This procedure was repeated at each temperature below Ty (~ 110K for devices under 
test; see figure [2b inset), down to ~ 75 K. At T < 75 K, the high values of Vsw necessary 
led to irreversible alteration of the electrode geometry, resulting in asymmetric, distorted 
Vsw histograms. Near 80-90 K, I-V cycles with multiple Vsw events (figure [lb) were 
observed more frequently. Thus, the total number of switching events observed varied 
with T, even with a fixed number of I-V cycles at each temperature. To compare Vsw 
distributions at different temperatures, the distributions were normalized, plotted as 
of counts)/(max # of counts) vs Vsw ; where "max # of counts" is the number of events 
at the most probable Vsw and of counts" is the number of events at a certain Vsw- 
Figure [2|i is an example of normalized Vsw distributions in the 77 K-105 K temperature 
range. The measured widths of the Vsw distributions are not limited by temperature 
stability. 

As has been discussed elsewhere [H], the use of pulses is essential to minimize the 
role of self-heating once the system has been driven into the more conducting state. This 
self- heating and the short timescale[ll] required to raise the local temperature in the 
channel significantly makes it extremely challenging to determine directly whether the 
initial breakdown takesplace through the formation of a conducting filament or through a 
uniform switching; once a highly conducting path is formed, the whole channel rapidly 
becomes conducting through self-heating. The filamentary picture is certainly likely, 
based on other breakdown phenomena in solids, and the statistical variation in Vsw is 
consistent with the idea of a process involving run-to-run variability associated with 
local details rather than global material properties, but this is not definitive. 

The Vsw distribution at each temperature is characterized by two main parameters: 
the mean switching value Vsw = (X)£Li Vsw,i)/-^; where N is the total number of 
switching ev ents; and the width of distr ibution, calculated as a standard deviation: 
c(Vsw) = \J (J2iLi(V S w,i — V sw ) 2 )/N — 1. As expected, the Vsw(^) has the same In- 
dependence (see figure [2] b) as Vsw(^) in single I-V experiments described in previous 
publications [SI El- More interesting is the cx(Vsw) temperature dependence, showing 
broadening of the Vsw distribution as the temperature decreases (figure [2] c) . We note 
a deviation from monotonous temperature dependence of <t(Vsw) at 100 K, observed in 
several devices tested. This is a temperature well below Ty = 110 K (see figure[2b inset), 
where several physical parameters (resistance, heat capacity and magnetoresistance) 
change abruptly. 

This increase in cr(Vsw) as temperature decreases is rather counter-intuitive. One 
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might expect "freezing" of temperature fluctuations and decrease in thermal noise the 
temperature decreases and, thus, narrowing of Vsw distributions. The field-driven 
breakdown of the insulating state is an example of the "escape-over-barrier" problem, 
addressed generally by A. Garg [19]. Below T v , the (temperature-dependent) effective 
free energy of the electronic system is at a global minimum value in the insulating state, 
while the external electric field modifies the free energy landscape, lowering the free 
energy of another local minimum corresponding to the more conducting state. As the 
external field is increased beyond some critical value, the minimum corresponding to 
the more conducting state becomes the global minimum. The nonequilibrium transition 
to the conducting state then corresponds to some process that crosses the free energy 
barrier between these minima. At a sufficiently large value of the external field, the free 
energy has only one minimum, corresponding to the conducting state. 

As Garg showed, one may consider thermal activation over the free energy barrier as 
well as the possibility of quantum escape. This free energy picture predicts a broadening 
of the transition driving force (Vsw here) distribution as the absolute value of the driving 
force increases, consistent with our observations. This free energy picture has proven 
useful in studying other nonequilibrium transitions, such as magnetization reversal in 
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Figure 2. (color online) (a) Normalized Vsw distributions at different 
temperatures (77 K - 105 K). (b) Temperature dependence of the mean 
switching voltage, Vsw- Inset shows zero-bias R vs T plot demonstrating 
Ty ~ 110 K.(c) Temperature dependence of Vsw distribution width, cr(Vsw) 
(black circles). 
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nanoparticles [20] and nanowires [2TJ Pinning due to local disorder is one way to 
find increasing distribution widths as T — > 0, as seen in investigations of field-driven 
magnetization reversal in nanowires [2T| 122] . Unfortunately, quantitative modeling 
in this framework requires several free parameters and is difficult without a detailed 
understanding of the underlying mechanism. 

Qualitatively similar phenomenology (distribution of switching thresholds that 
broadens as T is decreased) is also observed in the current-driven superconducting- 
normal transition in ultrathin nanowires |24l 125] . In this latter case as in ours, self- 
heating in the switched state is of critical importance, as is the temperature variation of 
the local thermal path. Again, quantitative modeling using this self-heating approach 
would require the introduction of multiple parameters that are difficult to constrain 
experimentally, as well as detailed thermal modeling of the nanoscale local effective 
temperature distribution, and is beyond the scope of this paper. 

We also examined the dependence of the switching distributions on applied out- 
of-plane magnetic field, H. Vsw distributions (3000 cycles each) were collected 
consecutively at 12 magnetic field values: T (first) — > 0.2 T — > 0.4 T — > 0.6 T — > 
0.8T^1T^2T^3T^4T^5T^6T^0 T(last). Figure E] a shows the 
resultant Vsw distributions at 80 K at several selected magnetic fields. As can be seen, 
magnetic field shifts the Vsw peak to higher V values and narrows the Vsw distributions. 
To reassure that the observed Vsw shift is not from irreversible changes in the device, 
a control experiment returning to H = T was performed after experiments in all 
non-zero magnetic fields. The Vsw distributions at H — T initially [H = T (first)] 
and in the end [H = T (last)] are identical (see figure [3] a), meaning that observed 
changes in Vsw distributions (shift and narrowing) are indeed caused by the applied 
magnetic field. Figure |3]b quantifies the dependence of Vsw an d cx(Vsw)- It is clear that 
both parameters saturate as H is increased beyond 1 T, i.e.; further increases of H up 
to 6 T have no significant effect. Magnetic field of the opposite polarity (not shown) 
has exactly the same effect on the position and the width of Vsw distributions. Note 
that the shape of the distribution, in particular its asymmetry about the peak value 
of Vswj evolves nontrivially with magnetic field, becoming more symmetric in the high 
field limit. 

We consider whether this magnetic field dependence of the switching originates 
with some dependence of the bulk resistance or contact resistances, as this would alter 
the electric field distribution in the channel. Since Vsw scales with the channel length, L 
0, HI] i.e., as does the resistance of the channel (R ~ L), one might expect an increase 
in Vsw could result from an increase in the device resistance with applied magnetic 
field. However, Fe 3 4 has a negative magnetoresistance (MR) [231 12H 127]. Figure [4] 
a shows an example of normalized resistance dependence, R/R(H = 0T), on the out- 
of-plane magnetic field at 80 K. The resistance remains effectively unchanged up to 
~0.4 T and then decreases as \H\ increases. Thus, when Vsw an d cr(Vsw) experience 
the predominance of their changes upon H application (H < IT, see figure [3]), the 
resistance of the device either stays constant or decreases. In the H range when R 
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Switching voltage, V sw (V) Magnetic field, fiH (T) 

Figure 3. (a) Examples of Vsw distributions at selected magnetic fields 
(T=80 K).(b) Magnetic field dependence of the mean switching value, Vsw 
(red squares), and the width of Vsw distributions, a (Vsw) (black circles). 
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Figure 4. Dependences of the resistance (R/R(H = 0T)) (a) and the 
magnetization (M/M s ) (b) on the out-of-plane magnetic field applied. 

experiences significant changes (see figure H] a), Vsw and cx(Vsw) remain essentially 
unchanged (figure [3]). Therefore, the shift of Vsw m the presence of H does not originate 
from the change in the resistance value of the device. 

Another Fe3C>4 film parameter effected by H is the magnetization of the film. 
Figure H] shows the normalized out-of-plane magnetization, M/M s , as a function of 
the out of plane H, where M is the magnetization of the film and M s is the saturated 
magnetization. This data is consistent with prior measurements on magnetite films |28] 
While we do not know the microscopic arrangement of M in the film in the absence of an 
external H, magnetostatic energy considerations mean that M under that condition lies 
in the plane of the film. The H range over which M is fully reoriented out of the plane 
(up to 1 T) matches the H range of changes in the position of Vsw and <t(Vsw, T) (fig. [3] 
b). This suggests (though does not prove) that the switching kinetics parameters Vsw 
and cr(VsW)^), and therefore the stability of the gapped, low temperature, insulating 
state is tied the magnetization direction of magnetite films. 

This observation is intriguing because it is not clear how the nonequilibrium 
breakdown of the low temperature state would be coupled to the magnetization. Possible 
factors include magnetoelastic effects such as magnetostriction [29] (~ parts in 10 4 
per Tesla) affecting the tunneling matrix element between B-site iron atoms; and 
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spin-orbit coupling playing a similar role[30j. There have been reports of significant 
magnetoelectric and multiferroic effects in magnetite [HIE], and a recent calculation [30] 
argues that these originate through the interplay of orbital ordering and on-site spin- 
orbit interactions of the B-site electrons. In this picture, reorientation of the spin distorts 
the partially filled minority-spin orbitals occupied on the B-site (formally) Fe 2+ ions. 
Such a distortion would be a natural explanation for the observed correlation between 
M and the kinetics of the electric field-driven breakdown of the ordered state, which 
directly involves the motions of those charge carriers. It is unclear how this kind of spin- 
orbit physics would explain the evolution of the Vsw distribution, however. It would 
also be worth considering whether there is any correlation between the characteristics 
of the switching distributions reported here, and the recently observed glassy relaxor 
ferroelectric relaxations in bulk magnetite crystals |3Tj. 

Additional, detailed experiments as a function of directionality of H, M, and 
crystallographic orientation should be able to test these alternatives. With the existing 
(100) films, studies of Vsw and cr(Vsw,T) as a function of H in the plane as well 
as perpendicular to the plane should be able to access the tensorial form of the H 
dependence. Comparison with appropriately directed M vs. H data as a function of 
temperature would be a clear test of whether the observed agreement between if -field 
scales (in Vsw an d reorientation of M is coincidental. Further measurements on films 
grown with different crystallographic orientations would serve as a cross-check. It is 
important to note, however, that the acquisition of such data is very time intensive 
due to the need to acquire many thousands of switching events. In turn, there is a 
companion requirement of extremely good device stability, to avoid irreversible changes 
in the metal configuration over the thousands of switching cycles. 

4. Conclusions 

We have studied the statistical distribution of the electric field needed for breakdown of 
the low temperature state of FeaO^ The distribution of critical switching voltages moves 
to higher voltages and broadens, as T is reduced. This broadening is consistent with 
phenomenology in other nonequilibrium experimental systems incorporating disorder 
and thermal runaway effects. The breakdown distributions are altered by modest 
external magnetic fields normal to the film, suggesting a need for further experiments 
to understand the connection between magnetization and breakdown of the correlated 
state. 
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